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ABSTRACT 


A method for solving a general class of three-dimensional 
boundary layer flows is developed. In the development, Levy- 
Lees variables are extended to three-dimensions and equations 
are placed in these similarity variables. An implicit finite 
difference scheme which is stable for negative transverse 
velocities is used to solve these equations. The method 
developed is applied to obtain solutions for sharp and 
spherically blunted circular cones at angle of attack. 
Longitudinal and transverse distributions are presented for 
these cases. Good agreement is found with the results 
obtained by other numerical schemes and the experimental data 
of Tracy, for sharp circular cones at angle of attack. For 
spherically blunted cones at angle of attack, the results are 
in good agreement with axisymmetric sphere results up to the 
region where spherical symmetry holds. 
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I. 


INTRODUCTION 


The importance of three-dimensional effects in boundary 
layer theory has long been recognized, as indicated by the 
reviews of the literature by Cooke and Hall [1] , Mager [2] 
and Stewartson [3], Due to the complexities of these flows, 
the majority of the literature has been concerned with 
approximate methods or with treating limiting cases where the 
boundary layer equations reduce to quasi- two-dimensional 
form. Yen and Thy son [4] have used integral methods and have 
shown reasonable agreement with experimental heat transfer 
data. By far, the most popular method of solution has been 
developed through use of small cross flow theory, which reduces 
the equations to a quasi- two-dimensional form. Chan [5] has 
shown remarkable agreement of the small cross flow theory with 
experimental data. 

As compared to two-dimensional boundary layers, relatively 
little work has been done on the computation of generalized 
three-dimensional boundary layers even for the simple case of 
the flow over axisymmetric bodies. Most of the work done in 
this area is restricted to only the windward plane of symmetry, 
e.g. Libby [6] and Watkins and Blottner [7], Another case 
which has received considerable attention is that of the 
three-dimensional boundary layer on a conical body. The 
general approach here is to transform the governing equations 
so as to eliminate dependence on the streamwise independent 
variable as has been done by McGowan and Davis [8] , Dwyer [9] 
and Boericke [10] . 
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Dor f.l.J.1 has used an explicit method developed by 
Haotx | 12] for solving the full three-dimensional boundary 
layer equations through use of zones of influence and 
dependence and the method has been applied to solve spherically 
blunted cone problems . 

The approach used here will be very similar to that of 
McGowan and Davis [8] in as much as the finite difference 
scheme and the general solution technique is concerned. 

However, here we will use Levy Lees 'type variables rather 
than working with the Grocco type variables of ref. [8]. 

An implicit finite difference scheme suggested by Krause [13] 
will be used so that regions of negative transverse velocities 
can be handled with the stability problem practically 
eliminated. 

The general three-dimensional boundary layer equations 
will first be specialized to Geodesic co-ordinates. A 
similarity type of transformation will then be introduced and 
the equations will be put in Levy Lees type of variables. 
Introduction of these variables will result in the correct 
similarity form of the equations, dependent only upon the 
normal transformed co-ordinate at stagnation points and sharp 
edges. Since the dependence on one co-ordinate drops out 
along the lines of symmetry, we can generate two planes of 
starting solutions originating at the sharp edge (or the 
stagnation point) from which the entire solution can be 
obtained using step-by-step integration. 
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The purpose of this investigation is to develop the 
three dimensional analogue of the Levy Lees variables that 
can be used for a general class of three-dimensional boundary 
layers and apply these to study the flow around sharp and 
spherically blunted circular cones at angle of attack. 
Calculations will first be made for 10° half-angle sharp 
circular cone at an angle of attack and the results will be 
compared with numerical calculations of McGowan and Davis [8] . 
The heat transfer data will be compared with experimental 
data of Tracy [14] . Flows around a spherically blunted cone 
of 10° half angle and at 4° and 8° angle of attack will be 
investigated and the results will be compared with the 
numerical results of Popinski * [15] . The inviscid data for a 
sharp circular cone has been taken from numerical calculations 
of Jones [16] while for blunted cones the data has been taken 
from numerical calculations of Rakich and Cleary [17] . 
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II. FORMULATION OF THE PROBLEM 


The problem will now be put in a form suitable for finite- 
difference analysis. The governing equations will first be 
written in generalized body-oriented co-ordinates and then 
specialized to the Geodesic co-ordinate system. Finally, a 
Levy-Lees type of transformation will be introduced to put the 
equations in the similarity variables so that similar solutions 
can be obtained at proper points. As a result of this trans- 
formation, the boundary layer growth is practically eliminated. 

The governing equations for three-dimensional compressible 
boundary layer flow in generalized body-oriented co-ordinates 
as shown in Figure 2.1, are given by Mager [2] as: 

5 -momentum 


_ lii + Y_ m + w_ + H v 
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n -momentum 


u 


h. 


av + v_ 


av 


w 


35 h 2 3n h. 


av 


ah 


-2 

u 


h K 3n 

n l 2 


, aH~ 

1 u v 2 
h l h 2 


l ap l a r p av, 

an “ -v- ac it as J ' 

ph 2 ph 3 h 3 


( 2 . 02 ) 


- 4 - 



Energy 


u 



at v at 

9 5 tt 9n 

2 


w at 




[ 


u ap 

V 5 


i£i 

9n J 


S Bh 3 2 


[(~) 2 
1 v a? ; 


(—) 2 ] 

v ac ; J 


+ 


l 




a_ 

3C 




t 


and Continuity 


(2.03) 


a 

n 


(h 2 h 3 pu) 


an {h i h 3 V) 


k (h l E 2P W) - 


(2.04) 


In addition, the equation of state and the viscosity law may 
be written as: 

Equation of State 


p = p R T 


(2.05) 


and Viscosity Law 


p = y(T) (2.06) 

Equations (2.01) through (2.06) are now nondimensionalized 
according to the following scheme: 
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( 2 . 07 ) 


( 2 . 08 ) 


/Re , 


( 2 . 09 ) 


( 2 . 10 ) 


( 2 . 11 ) 


/Re , 


( 2 . 12 ) 


( 2 . 13 ) 


( 2 . 14 ) 


-2 /K 
u /R 


( 2 . 15 ) 


-6 
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(2*16) 
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CO oo 


2 ' 


and 



(2.16a) 


where 



(2.17) 


In this scheme, Reynolds number is used in the nondimen- 
sionalization of h^ and w purposely so that any explicit 
dependence of the final equations on the Reynolds number is 
eliminated. 

Introducing equations (2.07) through (2.17) into the 
equations (2.01) through (2.06), we get the following set of 
nondimensional equations: 

£ -momentum 
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v 


h l h 2 



uv 

h l h 2 



1 8p 1 3 , y 9_Uv 

ph 1 a e , P h 3 ac K h 3 ac ; ' 


(2.18) 


n -momentum 
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Energy 


h x 35 h 2 3 n h 3 9 c YP ^ 95 h 2 9n ; 
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Continuity 


|j (h 2 h 3 pu) + |^- (h 1 h 3 pv) + |^- (h 1 h 2 pw) = 0 


State 
P = Pt , 
and 

Viscosity Law 
y = u (t) . 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


In general, an infinitesimal arc length is given by 

d£' 2 = (h-jdO 2 + (h 2 dn) 2 + (h 3 dc) 2 (2.24a) 

If we take d5 to be the infinitesimal arc length along the surface in 
the streamwise direction, h^ = 1 and assume the classical thin boundary 
layer approximation, h 3 = 1 equation (2.24a) may be written as 

d£' 2 = (d?) 2 + h 2 2 (dn) 2 (dc) 2 (2.24b) 
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and we let 


h 2 = h 


(2.25) 


Using (2.24) and (2.25), the governing equations (2.18) 
through (2.21) take the following form, whereas the equation 
of state and viscosity law are still represented by equations 
(2.22) and (2.23) respectively. 

5 -momentum 
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and 


Continuity 
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Introducing the following transformation in the normal 
direction 
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Let us define a velocity-like term V as 


V = - [— /2j I FdZ] 

q 3 £ Y 

e o 


- {M. |_ [_SL {M. f GdZ] 
q 3 n w h 

e o 


+ h p W 
K w w 


/2± 


(2 


where 


q = p p W h 
M w e. e e 


(2 


Then the transformed continuity equation is 
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and the normal velocity component is given by 
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Using this transformed continuity equation, the convective 
operator takes the form 


u L+VL + w L = uF^ + |- + v q fy 
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From similarity considerations, <[>(5,n) is found to be of the 
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!.40) 

form 
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(2.41) 


. = / p e P e W e h + S(n) 

o 

and W is chosen to be equal to u in our case. 

6 6 

Using the transformation of the normal co-ordinate (2.31) 
and the convective operator from (2.40) and letting 


£ = 


P P 
e e 


the governing equations take the form: 
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and Continuity Equation 
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The coefficients A^-A-^ depend only on the inviscid quantities 
and geometry and are defined as follows: 
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For the starting point where £= 0, the values of these coefficients 
are determined by a limiting process. 

The boundary conditions on F, G and 0 are at Z = 0: 


and 



(2 . 60a-c) 


where t is the wall temperature which is assumed specified and 
as Z -*■ 


F = 1, 

G = 1, 

and 

9=1. (2.61a-c) 
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The viscosity coefficient will be given by Sutherland's law: 



(t r +C) t* 3 / 2 
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where C is the Sutherland constant which is equal to 198. 6°R 


for air. 
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Nondimensionalizing the temperature by t r = we get 
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From (2.65), we get 


1 d«. _ (C*-6) 


i do 20 ( 9+C* ) * 


(2.67) 


The nondimensional skin friction coefficients can be shown 
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and 


, 9G p e v e h 
C = 2u V % ry 

f n e e 82 m 


(2.69) 


and the nondimensional heat transfer rate can be shown to be 
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Finally, the physical normal distance can be expressed as 


Z 
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e e 
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This formulation is very similar to that of McGowan and 
Davis [8] except for the fact that similarity variables are used 
here . 
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III. FORMULATION OF THE PROBLEM IN 


FINITE DIFFERENCE FORM 

The numerical scheme used to solve the coupled nonlinear 
second-order partial differential equations (2.42) through (2.45) 
is presented here. In principle, these three equations are 
replaced by a set of consistent linearized algebraic equations. 
The linearization is achieved by updating the coefficients 
appearing in these equations through iteration. This set of 
equations is of tridiagonal form which can be solved easily for 
boundary conditions of the type (2.60) and (2.61). 

The finite difference scheme used here is the one developed 
by Krause [13] and it can be Used for negative transverse 
velocity regions. This being an implicit scheme, problems of 
stability and mesh size are greatly reduced. The set of finite 
difference equations so developed is uniformly valid to second- 
order in the mesh spacings of the finite difference grid. The 
general grid system used is shown in Fig. 3.1. The equations 
are written at the mid point (C,n) and solved for the quantities 
F, G and 0 at the point (2,n) , assuming that these quantities 
are known from previous calculations at the points (l,n) , (3,n) 

and (4,n). After solving the energy, ^-momentum and n-momentum 
equations, the continuity equation, which is a first-order 
partial differential equation, is solved at the mid point (C, n ) 
subject to the boundary condition 

at Z = 0 V (C,n ) = 0 . 
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An iterative scheme is used on all the 4 governing 
equations, until the value of skin friction differs very 
slightly from the previous calculations. 

Three separate cases will be handled; a case in which 
similarity exists in both surface variables £ and n , a case 
in which similarity exists in one of the surface variables, 
say £ or n and a scheme to be used in general case. Before 
we discuss these cases, it will be more convenient to 
rearrange the equations (2.42) through (2.44) as follows: 
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Thus it is clear that if W represents F, G or 0> whichever 
is applicable, then all of these three equations can be 
written in the general form 
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where the coefficients A^ n , A 2 n , A 3n , A^ n and A 5 n are given by 
For the g -momentum equation 
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A ..V 

5n l 


(3.05) 


(3.06) 


(3.07) 


(3.08) 


(3.09) 


For the p-momentum equation 


1 d?. 3 8 V 

A in i de az i ' 


(3.10) 
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f 


(3.11) 


A 


(A2G+AqF+A^qF ) 


2n 


A 


^ A 2 + A 8 + A 10 ^ 


3n 


A. 

4n 

and 


„ F 
A 7 Z ' 


A 5n = " A 


4 Z * 


and for the energy equation 


A _ 1 dn 98 _ Pr „ 

In Z d0 9 Z z ' 


fl 2„ - - T- A 11 (F - G) ’ 


A 


3n 


Pr [A 13 ( a7 -) + A 14 ( gz ) ] 


4n 


Pr A 7 T ' 


and 


A 


5n 


pr A 4 ! 


(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 


The following discussion will refer to the general form 
(3.04) of the governing equations, rather than the original 
equations . 
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Case I . Solution at 5=0 and n=0 

This refers to the case where similarity exists in both 
the 5 and n directions. So the governing equations reduce to 
ordinary differential equations. In this case (l,n) coincides 
with (2,n) and (3,n) coincides with ( 4 , n) because of the 
similarity in n-direction whereas ( 2 ,n) coincides with ( 3 ,n) 
because of similarity in 5 -direction. So the equations are 
written at the point ( 2 ,n) and the solutions are also 
generated at this point only. The non- zero partial derivatives 
in equation (3.04), which are now actually total derivatives, 
are of the form: 


ItL ° M 2 ' n 2 AZ ~ WZ ' n ~ 1 + °<* Z > 2 ' 

2 , n 


(3.20) 


and 


3 2 W 


3Z 


2 , n 


W 0 - 2W 0 + W~ , 0 

■= 2 L n+1 — 2/ - £ 2 L n ~ 1 + 0(AZ) 2 . 

AZ 2 


(3.21) 


In this case, the governing equations are of the form 


3 2 W 


3 Z 


+ A 


3W 
In 3 Z 


n 


n 


+ A 0 W + -A, =0 

2n n 3n 


(3.22) 


2 

Substituting for — — ^ and from (3 . 20) - (3 . 21) into (3.22) v/e 

3Z 2 82 

get the tridiagonal form. 


A n W 2 ,n -1 + B n W 2 , n + C n w 2 , n+i = D n 


(3.23) 
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where 


A - 
n 



A 


— J 

2 A Z 


/ 


B 

n 



+ A 2„> 


/ 


c 

n 


<-V ^s> 

AZ 2AZ 


/ 


and 


D = 
n 


A 


3n 


(3.24) 

(3.25) 

(3.26) 

(3.27) 


Case II . Solution at g=0 or n=0 

This refers to the case where similarity exists in either 
the ^-direction or the n-direction. A Crank-Nicolson finite 
difference scheme is used for the case C=0 , while a truely 
implicit-finite-difference scheme is used for n=0. 


a. Similarity in the g-direction 

In this scheme (Figure 3.2a), the equations are written 
at the point (C,n) and solved for the quantities F, G and 0 at 
the point (2,n) knowing the solution at the point (l,n). The 
partial derivatives can be written in the form 


W — XW~ + (1— x) W. + 0 (An) t 

C, n 2 , n 1, n 


(3.28) 


9W 

az 


X(W 2,n + l - + - W l,n-1> + 0(4 "> ' 


C , n 


2 AZ 


(3.29) 
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- 2 W 
d W i 

zp' 

r ) L 


x[W 2,n+i 2W 2,n +W 2,n-l ] + (1 X)[W l,n+l 2W 1 ,n +W l ,n-l ] 

,2 


C , n 


AZ 


+ 0 (An ) § 


(3.30) 


and 


9 W 
9n 


C, n 


w 0 _ - w, 0 

-ki kil l + o (An ) 2 

An 


(3.31) 


where X - Crank-Nicolson number. 

Substituting these in the equation (3.04) and rearranging it 
in tridiagonal form, we get 


A_ = (-4- ii2)x , 


n 


AZ 


2 2AZ 


(3.32) 


B - ( + A 9 ) A + / 

n aZ 2 2 , n An 


(3.33) 


c = (j_ + ^12) x 
n '2 2AZ ;A ' 
A L 


(3.34) 


and 


‘ D n - A 3,n + ^ <777 + Til* W l,n + 1 

AZ 


+ (l-x) (- + A~ )W 


AZ 


2n 1 , n 


+ (1-X) (- 


In 


) W 


5 n 


W, 


AZ 


2 2 a Z 7 l, n -l An 1 , n 


(3.35) 


b. Similarity in n~direction 

Here the equations are written at the point (C,n) and 
solved for values at the point (2,n), knowing the solutions at 
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the point (3,n). The finite-difference grid is shown in 
Figure 3.2b. The partial derivatives become 


V7 


W~ + W 0 0 

= 2 - n 3 ' n + 0(4 5 ) 2 , 

^ t n ^ 


(3.36) 


9W 
9 Z 


(W 


2^n+l 2 / . n . Vi A r> l^n+l 3 , 'n-lj , + 0(AC) 2 , (3.37) 


C,n 


- W~ , + W- 
2_f_r 

4 AZ 


9 2 W 


9 Z‘ 


C,n 


2AZ 2 


+ 0(AC) , 


(3.38) 


and 


9W , W2,n W3,n + 0 (AC ) 2 • 

34 C ,n 


(3.39) 


Substituting these in the equation (3.04) and putting it in 
the tridiagonal form, we get 


A 

n 



In 

4AZ' 


/ 


B 

n 





C 

n 



In 
4 AZ 


t 


(3.40) 


(3.41) 


(3.42) 


and 
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A 


" n n A J,n +W 3,n+1 ( J^2 + TEZ 


- 1 * + W 3,n <- ~5 


1 A 2n 
' ^ 


A 


4n, 


AZ 


AT, 


A 


+ W 3,n-1 ( 


in, 


2 A Z' 


4 AZ J 


(3.43) 


Case III . General Case, £70 , n70 

Here the equations are written at the point (C,n) and 
solved for the values at the point (2, n) knowing the solution 
at the points (l,n), (3,n) and (4,n). The finite difference 

scheme is sketched in Figure 3.1. The partial derivatives can 
be written as 


W 


WL + W, 0 

= , - 3 ^ . + 0 (A 2 ) , 

c,n 2 


aw 

az 


aw 

az 


2 ,n 


aw 

az 


+ o (a 2 ) 


C,n 


(3.44) 


VI - w + w - w 

_ V 2 ,n +1 W 2,n-1 + W 3,n+1 W 3,n-1 , a/a 2, 

— ?AZ + 0(A ) 


(3.45) 


a 2 w 


2 

aw 


az‘ 


az 


2 ,n 


a 2 w 

az 2 


+ o(a 2 ) 


! C,n 


W- , -2W~ +W, _ ,+W, -2W a +W,. . 0 

= 2 ; n+l 2 f n 2 f n 1 3 , n+X 3 f n 3 r n 1 4 . q ^ 2 j (3 46) 

2AZ 2 


aw 


C,n 


w o « “ w , „ 0 

-^d? Id?. + 0 ( A 2 ) 

A? 


(3.47) 


and 
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+ 0(A 2 ) 


(3.48) 


'> VI i 

') n I 


(W 


2 , n 


-w, n ) 

I / n 


+ (W 


4 , n W 3 , rJ 


C, n 


2 A x] 


Substituting these into the governing equations and rearranging 
in the tridiagonal form, we get 


A 


A_ = [- 


n 


B 


2aZ' 

A 2n 


In 
4 aZ 


] , 


n 


[■ 


1 , A 4n A 5n-, 

+ — — + ] 


AZ 


AC 2a 


n 


(3.49) 


(3.50) 


C = 
n 


2 A Z‘ 


+ _iHi 
4aZ ' 


(3.51) 


and 


D " = A3n+< II? + ^ )W 3,n + l 


+ (- 


2 AZ 


In 


)W- 


2 4AZ' 3 , n-1 


5n 

2 A n ' W l , n 


+ (- 


2n 


AZ‘ 


A 4n A 5n 
AC 


) W 


2hr\' 3 , n * 


(3.52) 


The solution of the equations of the tridiagonal form 
(3.23) can be obtained by the standard methods for solving a 
linear tridiagonal algebraic set of equations. Knowing the 
values of the dependent variables F, G and 0 at n=l and at 
n=iE (corresponding to the wall and the outer edge of the 
boundary layer respectively) from the boundary conditions (2.60)- 
(2.61), each of the sets of equations represented by the general 
form (3.23) can be solved for IE-2, evenly spaced points from 
n ~2 , one step away from the wall, to the point n=IE-l, one step 
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away from the outer edge of the boundary layer by applying 


W 


2 ,n 


E 

n 


W, 


n+1 


+ F n 


(3.53) 


which is the general solution of (3.23) as justified by 
Richtmeyer [18] . 

Making use of the boundary conditions (2.60)-(2.61) and 
the fact that (3.53) should hold independent of the step size, 
we get 

El = 0 (3-54) 


for F, G and 6 and 
FI = 0, 0, t w 


(3.55) 


for F, G and e respectively. 

It can be shown that the following recursion relations 

hold 


-C 


E n = 


n 


A E - + B 
n n-1 n 


F = 
n 


D - A F , 
n n n-1 

A_E , + B 
n n-1 n 


(3.56) 


(3.57) 


Thus, E n and F n can be calculated knowing El and FI in 
addition to the coefficients A^ , B n , C n and D n , and using the 
above recursion relations. Making use of the outer edge 
boundary condition, W 2 n can be calculated at all the inner 
points of the mesh from the equation (3.53). 
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After solving the equations (3.04) in this manner, the 

coefficients A, , A 0 , A., , A„ , and A c appearing in the 
In ^n -> n 4 n 5 n 

governing equations are recomputed, updating the variables 
wherever needed and the equations are solved again. This 
process is repeated until the desired accuracy is achieved. 

General Solution Technique 

The solution technique used here is similar to the one 
used by McGowan and Davis [8] , i.e. a step-by-step integration 
is performed. As shown in Figure 3.3, two planes of starting 
solution are needed for such an integration scheme, along with 
the boundary conditions at the wall and at the outer edge of 
the boundary layer. 

At the intersection of these starting planes, the solution 
must be independent of both 5 and n - co-ordinates . In the case 
of sharp cone, this occurs at the sharp point while for blunt - 
body it occurs at the stagnation point. The solution technique 
for this case has been described under Case I in this chapter. 
Once the corner solution is obtained, the starting solution can 
be generated in n-plane at 5=0, using the method discussed 
ear li e r in this chapter under Case II. Also, the starting 
solution is generated on a line of symmetry with respect to 
the 5 -axis or along the windward streamline in our case. 

With the solutions known in these two starting planes, we 
proceed along constant 5 lines around the body (in the in- 
direction) and step down along the body (in the 5-direction) , 
once the solution in a particular 5-plane is completed. 
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IV. SOLUTION TECHNIQUE 


In this chapter we v/ill discuss the treatment of the 
inviscid data and the details of the solution technique 
employed herein. 

A. Inviscid Data 

The inviscid data is treated in much the same manner as 
by McGowan and Davis [8] . Thus we take the pressure distribu- 
tion from independent numerical solutions to the inviscid 
equations and then solve for the velocity distribution, such 
that it satisfies the Bernoulli Equation, longitudinal and 
transverse momentum equations ' and the isentropic relation, 
applied at the surface. 

The compressible Bernoulli Equation applied at the surface 

is, 


h = - 
° P, 


(—*y) 
v 1 


u 


V 


( 4 . 01 ) 


The ^-momentum equation applied at the surface is. 


3u v 3u 

u — — + r- — — — 
e ^ h 3n 


v 


ah 

35 


1 9P e 


p e h 9C 


( 4 . 02 ) 


The n~niomentum equation applied at the surface is, 


+ Is. !Ie + Ve 3h 

U e H h 3 n hac 


3P. 


P e h 3n 


( 4 . 03 ) 
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and the isentropic relation, assuming that the entropy is 
constant is, 


k = 


P Y 
K e ' 


(4.04) 


Differentiating (4.01) with respect to 5 


i “ -- 9P 3U 

n = Y Y -1 v-Y p Y ® + u ! 

Y — 1 Y e 95 e 95 


9v. 


+ V 


e 95 


and using (4.04), we get 


1 9P e 
P e 3? 


9u 


3 v. 


= u 


+ v 


e 95 e 95 


(4.05) 


Combining (4.05) with the 5-momentum equation and rearranging 


3 U 


e 


9 n 


9 

9 5 


(h v e ) = 0 


(4.06) 


which is the zero-vorticity equation applied at the surface. 

Note that if we differentiate the Bernoulli Equation with 
respect to n and combine with the n-momentum equation we will 
obtain (4.06) again. Thus, if we solve the Bernoulli Equation, 
isentropic relation and zero-vorticity equation, the solution 
v/ill automatically satisfy the 5-momentum and ri _ niomentum 
equations . 

For the special case of a sharp cone, equation (4.06) 
reduces to the simple form 
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(4.07) 


<) u 

o 

V ” V'- 

(i 'J 0 

because for this case 
h = C ■ 

8v e 

and — — = 0 since there is no variation in stream properties 
along rays, for conical flow. 

The general technique of generating the inviscid data from 
a specified pressure distribution has been described in detail 
by McGowan and Davis [8] . The same procedure is followed here 
for the sharp circular cone at an angle of attack and the 
pressure distribution for this case is taken from numerical 
calculations of Jones [16] . 

A slightly different approach is used for the case of a 
spherically blunted cone due to the nature of the problem* The 
pressure data for this case is taken from the numerical 
calculations of Bakich [17] . Up to the point where spherical 
symmetry exists, the data is specified in wind co-ordinates 
centered at the stagnation point (Figure 4.1) and from there 
onwards, the body-oriented co-ordinates, centered at the 
frontal point on the sphere are used for specifying the data. 

The data is specified in 7 planes evenly spaced around the 
body. 

For the symmetric sphere, the data is fitted with least- 
squares method to an even 14th order polynomial of the form 

2 4 14 

P = P Q (1 - P 2 S - P 4 S P 14 S ) (4.08) 

> 
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where S is measured from the stagnation point. 

in order to perform the 3-D calculations in body-oriented 
co-ordinates, it is necessary to find an expression for the 
pressure distribution as a function of S' measured from the 
frontal point O' as the origin. 

An 11th order polynomial of the form 

P = P 1 + P 2 S' + P 3 S' 2 + + P 11 S' 10 (4.09) 


is found for each of the seven planes in which the pressure data 

is available. For increased accuracy, separate polynomials are 

found for the three regions marked I, II and III in Figure 4.2, 

along the body. The pressure derivative in the longitudinal 
3P 


direction 


9C 


, is found by differentiating the expressions 


(4.08) or (4.09) analytically, in their respective regions. 

9P e 

Knowing the pressure distribution (P and ) , the 

O C, 

procedure for generating the inviscid data is the same as for 
the sharp circular cone, except for the fact that we use a 
more general equation (4.06), rather than the equation (4.07) 
(which is applicable only for conical flows) , in addition to 
the Bernoulli Equation and the Isentropic relation. 

Once the inviscid quantities are calculated at a convenient 
number of points around the body, a five-point Lagrangian 
interpolation is used to find the desired quantity,, if the point 
of interest lies, between the planes where inviscid data has been 
generated. 
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B. Solution to the Problem 

As described under the section "General Solution Technique , 
the planes of starting solutions can be generated, making .use 
of the fact that a similar solution exists at the windward 
streamline and then step-by-step integration can be performed 
making use of the boundary conditions at the wall and at the 
outer edge of the boundary layer. For the sharp cone, this 
procedure does not present any difficulties, but for blunted 
cones, some special attention is required to apply the 
technique. 

As mentioned earlier, the starting solution is obtained 
at the stagnation point in the wind co-ordinate system 
(Figure 4.1). The symmetric sphere calculations are carried 
out up to S=2ct, in a two-dimensional manner. The quantities 
F, G and 0 obtained from symmetric sphere calculations are 
stored as double arrays so that we can find their values as a 
function of the co-ordinates 5 and Z. 

Once we reach S=2a, we switch over to the body-co-ordinates. 
The starting solution in the plane X' — 1 — cos a, of constant 
5 around the body in this new co-ordinate system is then 
determined by interpolating the results of the symmetric sphere 
calculations. Note that in this new co-ordinate system we have. 


F' = F , 
G' = F , 

and 

0 ’ = 0 , 


(4.10a) 

(4.10b) 

(4.10c) 
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where ' refers to the quantities in the new co-ordinates. 

Once the starting solution in the plane £= constant is 
known , the general technique of step-by-step integration, 
described earlier is used to march along the body. 

A mention must be made of the manner in which the 
similarity parameter 4 . is handled in the new co-ordinates. 

The parameter <f> is found to be of the form 

K 

$ = f P e P e u e h 2 d£ . (4.11) 

o 

Let the calculation in a particular co-ordinate system be 
performed up to £ = £ . Then the expression (4.11) can be used 
to calculate <j> up to the point where £ =■ £ . 

Once we switch to a new’ co-ordinate system, the value of 
<j) should be calculated using the expression 

£ ' .. 

4> = 4> c + / p e n e u e h 2 d£ (4.12) 

where <j> q is the value of <f> in the starting plane £ = £ , 
obtained by using the equation (4.11), in the earlier co-ordinate 
system and integrating up to the point of interest in the 
starting plane. Knowing <fi , the values of <f> can be obtained in 
the new co-ordinates by direct application of the equation 
(4.12) . 

The value of h is determined using the expression for the 
arc length 
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dl' 2 = (dfj 2 + h 2 (dn) 2 + dZ 2 (4.13) 


and it is found that 


h = £ for sharp cone 


while h = local radius for the blunted cone. 
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V. RESULTS AND DISCUSSION 


Wo shall now discuss the results obtained for the 
particular cases solved. Comparison of heat transfer rates 
with experimental data will be made for sharp circular cones. 
The cases of sharp and blunted cones will be discussed 
separately. 

A. Sharp Circular Cone at Angle of Attack 

The case chosen here is one for which experimental data on 
heat transfer rate is available. The flow parameters for the 
case chosen are: 

M = 8.0, 

CO 9 

t. - S4.2°R, 

V*o = °' 46 ' 

R = 1718 ft^/sec^-°R, 

Pr = 0.738, 

C = 198. 6°R, 

e c = 10% 

a = 8°, 

and 

Y = 1.4. 

This particular case has been studied by Tracy [14] . The 
calculations are performed using 101 points evenly spaced in the 
normal direction and 31 points in the transverse direction, 
evenly spaced from n=0 to n=n 

max . 
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The inviscid pressure data of Jones [16] was used for these 
calculations. In Fig. 5.1 is shown the comparison of the 
pressure distribution used with Tracy's experimental data [14J 

7.95. Heat transfer rate distribution, nondimensionalized 
by q o = -0 . 0550 /ReC (which is the heat transfer rate for the 
circular cone at zero degree angle of attack) is shown in 
Figure 5.2. For comparison, numerical results of McGowan and 
Davis (8) and experimental results of Tracy [14J are plotted in 
the same figure. It is apparent that the results obtained are 
very close to the numerical results of McGowan and Davis [8], 
whereas agreement with experimental data is not as good, expecially 
near the windward and leeward regions. This could be due to 
differences in the pressure data of Jones [161 and Tracy [14] . 

in Figure 5.3, the longitudinal skin friction coefficient 
nondimensionalized by C f ^ = .40040 (where C f is the longi- 
tudinal skin friction coeHicient for the circular cone at zero 
d_ jree angle of attack) is shown. Numerical results of 
reference [81 are shown for comparison. The results obtained 
here are in good agreement for most of the region. 

Figure 5.4 represents the transverse skin friction 

coefficient nondimensionalized by C H^r-^ *». • 

y u f * Here also, the numerical 

results of reference [8] are presented ^ 

presented for comparison and they 

show a close agreement with the results obtained here. In 

conclusion we can say that for the sharp circular cone used here, ' 

we are in good agreement with the numerical results of reference" 

[81 and experimental results of reference [ 14 ], 
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B. Blunted Cone at Angle of Attack 


The test quantities used for this case are: 

M =8.0, 

CO 

t = 418 . 87 °R, 

oo 

t /t = 0.6, 
w o ' 

R = 1718 ft 2 /sec 2 =°R, 

Pr = 0.7, 

C = 198. 6°R, 

e c = 10°, 

a = 8° and 4°, 

and 

Y = 1.4. 

The inviscid pressure distribution is taken from the 

computer calculations supplied by Rakich [17] . This pressure 

2 

distribution, nondimensionalized by yP^M^, > is presented in 
Figures 5.6 (for a=8°) and 5.10 (for a=4°) for the windward and 
leeward planes. This inviscid pressure data, which is available 
for seven planes evenly distributed in the transverse direction 
from o=0 to T )=n max is represented by polynomial expressions for 
each plane and the pressure for intermediate values of trans- 
verse co-ordinate n are obtained by five point Lagrangian 
Interpolation. 

Figure 5.7 shows the heat transfer rate distribution along 
the body for the windward and leeward planes, for 8° angle of 
attack. Numerical results of Popinski [15] are also shown for 
comparison. Figure 5.8 shows the longitudinal skin friction 
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distribution for the same case and again the numerical results 
of Popinski [15] are shown. It is found that there is a good 
qualitative similarity between the two sets of results but 
numerical values are off by a few percentage points. This can 
be attributed to the manner in which the pressure data is 
treated in the two methods. It is noted that there is a slight 
kink m the pressure data near the sonic line from where 
onwards a different method of generating the data is employed. 
The kink results in a discontinuity in the pressure derivative 
8 p e /H and it has been found necessary to smooth the pressure 
distribution in order to avoid fluctuations in skin friction in 
this region. The final results will therefore depend on the 
degree of smoothing employed for the pressure distribution. 

Figure 5.9 shows the variation of heat transfer rate and 
the longitudinal and transverse skin friction distribution in 
the meridian direction for a station at S = 62° from the nose 
of the sphere. 

The results for a blunted cone at 4° angle of attack are 
presented m Figures 5.11 (heat transfer rate distribution), 

5.12 (longitudinal skin friction distribution) and 5.13 
(variation of heat transfer rate and longitudinal and trans- 
verse skin friction distribution in meridian planes) . 
Essentially, the same nature is exhibited by heat transfer rate 
and skin friction distribution as in the case of 8° angle of 
attack. But as expected, the difference in the values of heat 
transfer rate (or longitudinal skin friction) at the same 



longitudinal station for windward and leeward planes, in this 
case is smaller than in the case of 8° angle of attack. 
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dS^ = dr, 
df^ = hj dn 
dS^ = dc 

Fig . 2.1. BODY ORIENTED CO-ORDINATE SYSTEM 
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Z| 



Solution known 


Fig. 3.2a. FINITE DIFFERENCE SCHEME USED 
WHEN SIMILARITY EXISTS IN ^-DIRECTION 
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Fig. 3.2b. FINITE DIFFERENCE SCHEME USED WHEN 
SIMILARITY EXISTS IN INDIRECTION 
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Fig. 3.3. GENERAL SET UP FOR STEP-BY-STEP INTEGRATION 
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Fig. 4.1. CO-ORDINATE SYSTEM USED FOR SPECIFYING 
INVISCID DATA FOR SYMMETRIC SPHERE (x,y) 

AND FOR GENERAL 3-D BODY GEOMETRY (x',y') 
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Shock 



Fig. 4.2. SCHEMATIC DIAGRAM INDICATING THE 
THREE REGIONS WHERE SEPARATE POLYNOMIALS ARE 
USED FOR EACH REGION TO OBTAIN THE 
INVISCID PRESSURE DATA 
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F ig. 5.3. LONGITUDINAL SKIN FRICTION 
DISTRIBUTION ON SHARP CIRCULAR 


CONE AT 8° ANGLE OF ATTACK 


Present Results 


McGowan and Davis [ref. 8] 



n/n 


max 


Fig. 5.4. TRANSVERSE SKIN FRICTION 
DISTRIBUTION ON SHARP CIRCULAR 


CONE AT 8° ANGLE OF ATTACK 


Fig . 


5.5. SCHEMATIC DIAGRAM OF CO-ORDINATES 
USED TO PRESENT RESULTS FOR 
SPHERICALLY BLUNTED CONE AT AN 


ANGLE OF ATTACK 






Fig. 5.8. LONGITUDINAL SKIN FRICTION DISTRIBUTION FOR 
BLUNTED CONE AT 8° ANGLE OF ATTACK 
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F-ig. 5.9. VARIATION OF HEAT TRANSFER RATE, LONGITUDINAL 
AND TRANSVERSE SKIN FRICTION AROUND BLUNTED CONE 
AT 8° ANGLE OF ATTACK (S=62°) 
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Fig. . 5. 13. VARIATION OF HEAT TRANSFER RATE, LONGITUDINAL AND 
TRANSVERSE SKIN FRICTION AROUND BLUNTED CONE 


AT 4° ANGLE OF ATTACK (S=62°) 


